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Abstract. We present results of convective turbulent dynamo simulations including a coronal 
layer in a spherical wedge. We find an equatorward migration of the radial and azimuthal fields 
similar to the behavior of sunspots during the solar cycle. The migration of the field coexist with 
a spoke-like differential rotation and anti-solar (clockwise) meridional circulation. Even though 
the migration extends over the whole convection zone, the mechanism causing this is not yet 
fully understood. 
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1. Introduction 

The Sun shows an activity cycle with a period of eleven years, which is manifested 
by sunspots occurring on the solar disk. Every eleven years the number of sunspots has 
a maximum and around five years later a minimum. With each cycle the polarity of 
the global magnetic field changes sign and gives rise to the 22 years Hale cycle. The 
occurrence of sunspots does not only vary with time, but it also shows a strong lati- 
tudinal dependence. During solar minimum, the sunspots are more likely to emerge at 
higher latitudes and during the maximum at lower latitudes, close to the equator. This 
observed feature corresponds to an equatorward migration of the magnetic activity belt 
causing the sunspots. Producing the cyclic behavior in numerical simulations is chal- 
lenging. Either one neglects important backreaction of the magnetic field on the fluid 
motions and solves the problem in a kinem atic mean-field approach with widely varying 
parameterizations for the small scale s (e.g. Dikpati fc Charbonneau 19991 : Kapvla et al.l 
l2006t iKitchatinov fc Olemskovl 12012 ) or one uses direct numerical simulations (DNS) 
with orders of magnitudes too high viscosity and diffusivity, and too low stratification o f 
density and temperature (e.g. Brun et al. 2004 : Ghizaru et al. 20101 : Kapyla et a,l.ll20l"oh . 
The kinematic mean-field approach seems to be successful in reproducing observed fea- 
tures of the Sun. However, whether solar activity is generated by a distributed turbulent 
dynamo throughout the convection zone or in a thin layer below the convection zone is 
still very much under debate and cannot be answered just by mean-field modeling. It 
has turned out to be much more challenging to obtain solar-like magnetic a ctivity from 



DNS or large-eddy simulations: models eithe r show only weak m ean fields (jBrun et al 



20041). strong quasi-stead y large-scale fields ( Brown et al. 2010l). cycles but no migra- 



tion ( Ghizaru et aTl 2010l) . or cycle s and poleward migra tion ( Gilmanl 1983 : Kapvla et al 



2010HBrown et al.ll201ll) . Recently iKapvla et al.l (| 2 0121 ) were able to reproduce an equa- 



torward migration of the azimuthal and radial magnetic fields for the first time using 
direct numerical simulation of a convective dynamo. 
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Figure 1. Overview of the stratification. The normalized density p/po, pressure p/po, tem- 
perature T/To are plotted together with the specific entropy s/c p over the radius. The inset in 
logarithmic scale emphasizes the steep decrease of the pressure and density in the coronal layer. 



Another important observed feature of the Sun is its differential rotation. It is broadly 
believed that the solar differential rotation is an important ingredient in the solar dynamo. 
In the kinematic models the dif ferential rotation profile is u sed as an input parameter. 
Beside mean- field models (e.g. Kitchatinov fc Rudiger 19991 ). a solar- like differential ro- 
tation profile has not yet been self-consistently obtained in direct numerical simulations. 
Often the rot ation profile is so l ar-like in the sense that the equator rotates faster than 
the poles (e.g. Brun et al"1l2004l ; [Kapvla, et al.ll2011 ). but the angular velocity is constant 
on cylinders due to the Taylor-Proudman balance. Spoke-like rotation has not been seen 
unless a latitudinal temperature gradient is imposed at the base of the convection zone 
( Miesch et alll2006t) . 

In the following work, we present the results of a two-layer model, which combines a 
convection zone with a coronal layer. The model was originally developed to investigate 
coronal ejections generated by a dynamo undern eath the surface and the helicity flux 
due to the coronal envelope a nd the ejections, see lWarnecke fc Br andenbursd (j2010|) and 
Warnecke etafl (|201ll l2012allbh for details. In this work we focus on the magnetic and 
fluid properties in the convection zone. The dynamics of coronal ejections will be discussed 
in a followup paper. 



2. The model 

We use a similar setup as in IWarnecke et~al1 (|2012bl) . which we improve in two respects. 
Firstly, we apply more solar-like temperature and density profiles for the coronal layer. In 
the earlier work, the temperature in the corona was constant and identical with that in 
the photosphere. Now, we let the temperature rise above the surface reaching a constant 
value at a radius of r — 1.1 R, where R is the solar radius. The coronal temperature is 1.2 
times larger than that at the base of the convection zone. In Figure [1] the radial profiles 
of density, temperature, pressure, and entropy are shown from the saturated stage of 
the simulation. The entropy is negative in the convection zone with a negative slope. 
Close to the surface (r = R) the gradient changes sign resulting in positive values above 
r = 1.05 R. The density decreases exponentially in the corona (r > R). In the convection 
zone (0.7R < r < R) the stratification is initially isentropic. This results in a density 
difference po/p s = 14 in the convection zone and po/pt = 2.3 • 10 3 in the entire domain, 
where po — p(r — 0.7 R), p s = p(r = R) and pt = p(r = 1.6 R). The other improvement 
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is to use a turbulent heat conductivity \t in the entropy equation (|2.4[) to model the 
unresolved heat c ond uctivity. Similar subg rid-scale models have also been used by e.g. 



unresolved heat c ond uctivity, similar subg rid-scale models have also been used by e.g. 
Brun et al. ( 2004 ) and Kapyla, et al.l ( 20 1 lh to significantly increase energy transport due 



to the convective flux. We apply a constant value of Xt > which is similar to the molecular 
viscosity u, in the interval of 0.75i? ^ r ^ 0.97i?, and tends smoothly to zero above and 
below. 

We use a spherical wedge extending from r = 0.7 R to r — 1.6 R in radial direction, 
from 6 — 15° to 165° in colatitude, and from <f> = 0° to 90° in azimuthal direction. We 
solve the following equations of compressible magnetohydrodynamics: 

3 A 

* =ux B + nV 2 A, (2.1) 



dt 
Dlnp 



(2.2) 



Dt 

5^ = g - 2f2 x u+ - (J x B - Vp+ V • 2vpS) (2.3) 
Dt p 

T Sl = - V ■ ( KVT + XtpTVs) + 2^S 2 + M j2 _ r cool (r), (2.4) 
Dt p p 

where magnetic field is expressed by B = V x A, where A is the vector potential, 
and 77 is the magnetic diffusivity. Furthermore, u is the velocity, p is the density and 
g = —GMr/r 3 is the gravitational acceleration. The Coriolis and Lorentz forces are 
given by 2i~2 x u and J x B, respectively, where fi is the angular velocity and J = 
Pq x V x B is the current density. S is the traceless rate-of-strain tensor and v is the 
kinematic viscosity. The energy conservation equation is solved in terms of the specific 
entropy s. We include radiative transport in term of the diffusion approximation with 
heat conductivity K and model the unresolved turbulent transport with Xt- The fluid 
obeys the ideal gas law, p — (7 — l)pe, where 7 = c p /c v — 5/3 is the ratio of specific 
heats at constant pressure and constant volume, respectively, and e = c v T is the internal 
energy density, which defines the temperature T. The third and fourth terms on the 
rhs of Equation (|2.4j) describe viscous and Ohmic heating. The radial cooling function 
r c0 oi(f ) = To (eg — cl Q (r)/cl (r)) is defined such that it cools/heats the corona towards a 
predefined temperature profile c s o(r), which is plotted in Figure [TJ As an initial condition 
we use a stratification profile similar to that plotted in Figure Q] and a weak, random, 
Gaussian-distributed seed magnetic field in the convection zone. 

We apply periodic boundary conditions in the azimuthal direction. For the velocity 
field we use stress-free conditions on all other boundaries. For the magnetic field we use 
a perfect conductor boundary condition at the r — 0.7 R and both 6 boundaries; at the 
r = 1.6 R boundary be apply a vertical-field condition. The heat flux through the lower 
(r = 0.7 R) boundary is held constant and the temperature at the top is fixed. 

We employ the Pencil Code[J which uses sixth-order centered finite differences in 
space and a third-order Runge-Kutta scheme in time; see Mitra et al. ( 20091 ) for the 
extension to spherical coordinates. 



3. Results 

For this work we consider a single simulation, which is defined by the following prop- 
erties: the Coriolis number Co = 2il/u rms fcf is 14, the Reynolds number Re = u rias /uk{ 
is 20, and the magnetic Prandtl number Pr^/ = v/t) = 1. Here kf is the wavenumber of 

f http : //pencil-code . googlecode . com 
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Figure 2. Variation of the magnetic field in the convection zone. The azimuthal mean magnetic 
field and the radial mean magnetic field B r is shown at radius of r = 0.97 R (left panel) 
and r = 0.84 R (right). Dark blue shades represent negative and light yellow positive values. 
The dashed horizontal lines show the location of the equator at 6 = tt/2. The magnetic field is 
normalized by the equipartition value B cq . 



the energy carrying scales and is denned as kfR — 21, and u rms = y/3/2(v% + Ug) is the 
rms velocity at radius r — 0.84 R. 

Below the surface (r < R) the domain is convectively unstable. The resulting convective 
motions drive a turbulent dynamo and generate a large-scale magnetic field. The rms 
magnetic field B rms in the convection zone reaches 0.38 times the equipartition field B eq = 
yj Mo(p u 1 2 ms)} ''=o.84-R- In Figure [2] we plot the azimuthally averaged (in the following 
mean) radial and azimuthal magnetic fields as functions of time in terms of the turnover 
time t — ii rms fcf. Effectively r is an eddy turnover time, which varies between a few 
minutes and more than a month in the solar convection zone, depending on the depth. In 
the kinematic regime, the magnetic field pattern shows a cyclic poleward migration close 
to the surface (r = 0.97 R) for the azimuthal and radial fields B r at lower latitudes. 
The period is short, only a few tens t/r. In the non-linear regime, a new pattern is 
developing: equatorward migration starting at high latitudes overcomes the poleward 
migration at low latitudes. This equatorward migration shows a longer period of 100-200 
t/r. Additionally the pattern is not only visible near the surface, as it penetrates into 
the convection zone as shown in the left panel of Figure [2j The pattern does not look 
as regular as in iKapyla, et al. ( 2012 ). which can be explained by the interaction of the 
corona and convection zone. 

Initially the interior rotation profile is constant in cylinders with a positive gradient 
of angular velocity at low latitudes. In the hydrodynamic state, affected by rotation, the 
convective turbulence re-shapes the rotation profile within the convection zone, while 
in the convectively stable coronal layer only viscosity is at play. After a dynamically 
important magnetic field has been generated by the dynamo action, it influences the 
rotation as well. In the corona, where the density is much lower than in the convection 
zone, the Lorentz force can dominate over the Coriolis force, see Equation (|2.4[) . and have 
a greater influence on the angular velocity. In our simulation the velocities in the corona 
saturate not before t/r = 660. 

The differential rotation profile is solar-like, see Figure [3l The equator rotates faster 
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Figure 3. Meridional circulation and differential rotation zoomed in the northern hemisphere in 
the convection zone. Left panel: differential mean rotation profiles Q(r, 8)/Qo = u^j(r sin#) + 1 
overplotted with the poloidal velocity u p = u r + ug in white arrows. Right panel: the same as 
left, but without the poloidal velocity and emphasizing the mean differential rotation contours. 
The white-black dashed line indicates the surface (r = R). 



than the poles and the contours of constant rotation show a spoke-like pattern-especially 
at low latitudes. This is the first time a spoke-like profile ha s been found in DNS in 



combination with equatorward migration. iKapyla et al.l (|2012l ). whose setup is similar 



to ours in the convection zone, found that the contours of constant rotation are signif- 
icantly more cylindrical at similarly high Coriolis numbers. This could be due to the 
interaction between the coronal layer and the convection zone. There might be even an 
indication of a surface shear layer at latitudes 9 ^ tv/3. The meridional circulation shows 
a clockwise pattern-in contrast to the Sun. Instead of a poleward flow near the surface, 
as in the Sun, the flow is equatorward. This might be one possible explanation for the 
equatorward migration of the magnetic field. The field would then be dragged along with 
the fluid motion to the equator. But the magnetic field shows an equatorward migration 
throughout the convection zone, and the anti-solar meridional circulation seems to be 
shallow with a return flow in the middle of the convection zone. In our simulation the 
meridional circulation seems not to have a connection with the equatorward migration 
and the generation of the field. 

Another possible explanation of the equatorward migration can be the profile of the 
effective a which likely generates the mean magnetic field B. a has two contributions 



«k + a m = -\r c ijj ■ u + \r c j ■ b/p. (3.1) 



proportional to the kinetic helicity density u) ■ u and the magnetic helicity density j ■ b/~p, 
respectively, and where r c is a correlation time of the turbulence and can be expressed by 
the eddy turnover time r. In Figure[4j we show the two contribution as functions of radius 
r for two latitudes. In the convection zone a is dominated by the kinetic contribution 
and is influenced by a m only close to the surface. One reason behind the equatorward 
migration could be the different values of a at different latitudes. We observe a tendency 
for the a effect to have larger values at higher latitudes. The latitudinal profile of a 
could cause the mean field to migrate towards the equator. This has to be verified by 
comparison with simulations which show a poleward migration. 



Warnecke et al. 






Figure 4. Left panel: The two contributions to a plotted as functions of the radius r. Upper 
left panel: kinetic helicity density u> ■ u for two latitudes: 9 = 70° (black line) and 9 = 50° (red 
line). Lower left panel: current helicity density j ■ b/~p at 9 = 70° (black line) and 9 = 50° (red 
line). The values are normalized by Wj ms fcf. Right panel: a — — \tJ*) ■ u + ■ b/~p at 9 — 70° 
(black line) and 9 — 50° (red line). The dotted lines indicate the same as above, except in the 
southern hemisphere where it has the opposite sign, a are normalized by u rms . The dashed lines 
represent the surface (r = R) and the zero line. 



4. Conclusions 

In summary, we use a two-layer model that combines self-consistent convective dynamo 
action with a coronal layer, resulting in equatorward migration of the mean magnetic 
field. A solar-like rotation profile with s poke-like contours a t low latitudes is generated. 
This is different from recent results of iKapyla, et al.l (|2012l ) with a similar model, but 
without corona. In forthcoming work, we will study the mean magnetic field structure and 
differential rotation for different rotation rates and investigate the cause of equatorward 
migration in more detail. 
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